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ABSTRACT 

Context. Secondary eclipses are a powerful tool to measure directly the thermal emission from extrasolar planets, and to constrain 
their type and physical parameters. 

Aims. We started a project to obtain reliable broad-band measurements of the thermal emission of transiting exoplanets. 

Methods. Ground-based high-cadence near-infrared relative photometry was used to obtain sub-millimagnitude precision light curve 

of a secondary eclipse of WASP-4b - a 1.12 Mj hot Jupiter on a 1.34 day orbit around G7V star 

Results. The data show a clear >10o- detection of the planet's thermal emission at 2.2 fim. The calculated thermal emission corre- 
sponds to a fractional eclipse depth of 0.185^q q}3%, with a related brightness temperature in Ks of Tb = 1995 + 40 K, centered at 
Tc = 2455102.61 162+g;ggJJ^^ HJD. We could set a limit on the eccentricity of ^ cos a; = 0.0027 + 0.0018, compatible with a near- 
circular orbit. 

Conclusions. The calculated brightness temperature, as well as the specific models suggest a highly inefficient redistribution of heat 
from the day-side to the night-side of the planet, and a consequent emission mainly from the day-side. The high-cadence ground-based 
technique is capable of detecting the faint signal of the secondary eclipse of extrasolar planets, making it a valuable complement to 
space-based mid-IR observations. 
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1. Introduction Additional orbi tal and physical paramete rs w ere measured also 

tH- .11 1 by Gillon et al. (2009). Winn et al. (2009). and lSouthworth et alJ 

^ Photometry of extrasolar planetary transits allows to measure the ^ 2009). Recently, Beerer et al. ( 2011) performed photometry, us- 

Q planetary radii and density, and to infer the nature of the planet - • jj^^^ instrument on the Warm Spitzer, of the planet 

. gaseous or rocky The detection of the thermal emission during wASP-4b in the 3.6 and 4.5 ,,m bands. Their data suggests the 

g . the secondary eclipses when the planet passes behind the star wASP-4b atmosphere lacks a strong thermal inversion on the 

^ gives access to additional physical parameters like the brightness ^ .^^ ^^^^^^ unexpected result for an highly irradi- 

^ temperature, an important quantity tor the comparison with more ^^^^ atmosphere 

^ sophisticated planetary atmospherk models with chemistry and wASP-4b was one of the pilot targets for our program be- 

^> dynamics (i.e , Burrows et al. , 2QQ8). Whil e a few secondary ^^^^^ -^^ ^^^^ ^^^^ 1 ^j^^ ^-^^ ^ ^^^^ 

K> ■ e^^s es have been secure d from space (e.g | Charbonneau et alj transiting exoplanets, with a relatively large predicted secondary 

KN :l2QQilDemmgetal.' 2005), the detection of secondary eclipses ^^j^ utude of about a millimagnitude, placing it well 

. has been very challenging from the ground with only a ^j^j^j^^ ^^ ^^^^ detection. In Section 2 we discuss 

i t'ereralT Sg f^t":^^ *^ ^''^-^^^i-^ -^ly^i^- ^^^^i- ' P-^-^^ *^ 

I Kogers et al. | /UU^ | Lroli et al. || /UlUa||bt and reterences therein). ^^^^ secondary eclipse, and gives the measured brightness 

We obtamed sub-miUimagnitude precision ground-based re la- temperature for WASP4-b, comparing it with theoretical models. 

r.'^^P^ u"""^^"?' u °u : high-cadence near-infrared pj^^jj ^^^^^^ 4 the results and lists some conclu- 

(NIR) observ ations, which can det ect transits down to a few mil- ^[q^^ 
limagnitudes dCaceres et al .1120091: hereafter Paper I). In this pa- 
per we report the detection of a seco ndary eclipse of WA SP-4b. 

WASP-4b was discovered by Iwilson et aP (l2008l) . It is 2. Observations and data reduction 

a 1.12Mj hot Jupiter on a 1.34 day orbit around a G7V 
star, with a heavily irradiated atmosphere and large radius. 



The Ks band WASP-4b observations were obtained with 
the Infrared S pectrometer And Array Camera (ISAAC 



tne Inirared spectrometer And Array Camera (IbAAC 
iMoorwood etaD Tl998) at the ESQ VLT, on September 27, 



* Based on observations collected at the European Southern 
Observatory, Chile. Programme 083.C-0528. 2009. The observations were carried out in the FastPhot 

Correspondence to: cccdiCQrQ@diS\io.^\xc.c\ mode, as in Paper I. The InSb 1024xl024px Aladdin detec- 
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Fig. 1. Upper panel: The range of airmass covered during the 
run. Mid and bottom panels: The variation around the median 
values of the position of the centroid of the stellar psf on the de- 
tector. The positions of those images discarded because of point- 
ing jumps whose amplitude is larger than the plotted range are 
not shown, for the sake of clarity. 

tor was windowed down to 304x608 px, yielding a field of 
view of -45x90 arcsec. The field of view includes WASP- 
4 and one refere nce star - R. A 23 34 18.4, Dec. -4204 51.0 
(J2000). 2MASS (ISkrutskie etal.ll2006h lists for them respec- 
tively:/=11.2,if=10.8,and/^5 = 10.7magand 7=11.7,7^=11.3, 
and ^5 = 11.2 mag. 

This observation mode generates a series of data-cubes with 
a fixed number of frames. Each frame is an individual detector 
integration, and there is nearly zero "dead" time between the 
integrations. The cubes are separated by an interval of ~6 sec for 
fits header merging, file transfer, and saving on a hard disk. For 
these observations the cubes had 250 frames, corresponding to 
integrations of 0.6 sec each. In total, we collected 30,003 frames 
during a continuous 310min long run. 

Standard procedures for NIR data reduction as dark subtrac- 
tion and flat fielding were performed on the data set. We applied 
aperture photometry, using the IRAlQ package Daophot, to mea- 
sure the apparent fluxes of the two stars in the field. An aper- 
ture radius of 1.92arcsec was selected empirically to have the 
best compromise between the r.m.s. on the transit light curve, 
determined as ratio of the target- vs-reference fluxes, and the 
systematic noise introduced by a varying sky. The sky back- 
ground was estimated in a circular annulus with an inner radius 
of 3.10arcsec and outer radius of 7.5 arcsec, from the centroid 
of the point- spread function. 

Occasionally, during the observations the central peak of the 
stars exceeded the linearity limit of the detector of -6500 ADU. 
To avoid the non-linearity efl'ects that may distort the occultation 
depth measurement we discarded those frames. Moreover, at the 
end of the run the detector was moved yielding in jumps in the 
position of the stars. All the points after the first jump in posi- 
tion were also discarded. The starting points showed also a jump 
that could not be corrected, and these points were also removed 
to avoid including a systematic error to the determined param- 
eters. After cleaning the light curve, we leave a total of 24,070 
measurements, that were used during the subsequent analysis. 

The light curve shows a smooth trend during the period of 
observation, most likely due to airmass variation, important dur- 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 2. The eff'ect of rednoise on the obtained light curve is 
shown. The total noise for a given bin (solid line), is com- 
pared with the expected Poisson noise (dashed line). The resid- 
ual points for the complete light curve were used in this calcula- 
tion. 



ing the first half of the night, and due to target drift across the 
detector, which dominated the end of the run, as shown in Fig. 
[U To remove these eff'ects we modeled the base-line of the light 
curve with a polynomial that includes a linear dependency on 
time (0, airmass (sec z), and positions of the centroid of the star 
on the detector (xc and yc) as follows: 

fBline = + Cllt + a2 SCC Z + ^s^c + <34jc • (1) 

The coefficients in Eq. [T] were determined using the en- 
tire light curve, with the fitting procedure described below. The 
form of the polynomial was selected to flatten the out-of-eclipse 
lightcurve after multiple experiments. We also attempted to in- 
clude other parameters, i.e. the full-width-half-maximum of stars 
on each frame, but the correlation of the signal with this param- 
eter was negligible. 



3. Analysis and discussion 

Various parameters can be measured from the secondary eclipse 
light curve. The thermal planetary emission is proportional to the 
eclipse depth d, the orbital eccentricity e and argument of perias- 
tron w could be inferred through the analysis of both the eclipse 
central time Tc, and the secondary eclipse length Tocc- We create 
an occultation model following the algorithm for uniform illumi- 
nated transit li ght curve (i.e. neglectin g the limb-darkening con- 
tribution) from lMandel & Ag ol'('2002), and we scaled the model 
transit depth to fit the observed eclipse depth. The stellar and 
planetary parameters for WASP-4 and WASP-4b were adopted 
from Winn et al. (2009): period P=1.33823214d, planet-to-star 
radius ratio /?=0. 15375, semi-major axis a=5.473R^, orbital in- 
clination /=88.56deg, and stellar radius 7?*=0.9127?o. 

Past experience with NIR data showed that correlated noise 
could degrade photometric accuracy in these wavelengths. To 
evaluate the contribution of this kind of noise on our photometry, 
we analyzed the behavior of the light curve r.m.s. for diff'erent 
binning factors, after removing the best-fitting model. It shows a 
small deviation from the expected curve for pure Poisson noise 
(Fig. [2]) over range of binning corresponding to times of the order 
of tens of minutes - similar to the duration of the ingress-egress 
for the secondary eclipse. 

The best fitting model parameters were found wi th the mul- 
tidim ensional minimization algorithm AMOEBA (iPress et al.l 
Il992h . where the free parameters were the eclipse depth, cen- 
tral time, length, and the coefficients of the base-line in Eq. [T] 
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Fig. 3. The upper panel shows the 2-min bin uncorrected Hght 
curve, with the overlaid best fitting model considering the oc- 
cupation model and the base-line model. The corrected Ks light 
curve of WASP-4b is shown in the middle panel, including the 
complete data set, where the resulting rms in the out-of-eclipse 
points is 0.0069. In the bottom panel each data point corresponds 
to a 2-minute bin, along with its error, and the best-fitting oc- 
cupation model. The residuals for this curve are shown with an 
off'set of 0.006, and present a r.m.s. of 0.00072. A color version 
of this figure is available in the electronic form. 

on the 2-min binned light curve, yielding 131 data points. The 
function to minimize was the;^^ statistic, for a model that is: 



mod ■■ 



Bline 



(2) 



where Mocc corresponds to the occupation model described 
above. To determine the errors in the fitted parameters, we per- 
formed the bootstrapping procedure described in Paper I: we first 
subtracted the best fitting model from the data set, then we took 
the set of residuals and shifted the residual to become the i+V^ 
residual, and the last to become the first. Next, we added the re- 
ordered residuals to the best fitting model and run it through the 
minimization procedure instead of the data to obtain a new set 
of fitting results. This procedure was repeated until a full circle 
over the "good" points was completed, and the reported uncer- 
tainties were defined by the 68.3% level around the median val- 
ues of the distributions. Therefore, our errors properly account 
for the correlated noise in the data set. 

The final corrected light curve is shown in Fig. [3] Table [T] 
lists the individual flux measurements, with their time-stamps 
and Poisson uncertainties. The resulting parameters and their un- 
certainties for WASP-4b are listed in Table |2] 

The significance of our detection can be assessed from Fig.lH 
where the distribution of the relative flux (with respect to the 
reference star) during the eclipse and outside of the eclipse are 
compared. The bin width is equal to the calculated eclipse depth 
d = 0.185%, and each distribution curve was normalized for 
comparison purposes. A K-S test applied on these distributions 
shows they are identical with a ~ 99% of probability, but dis- 
placed exactly one bin. 
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Fig. 4. Distribution of the flux during the eclipse (dashed line) 
and out of the eclipse (solid line). The bin width is 0. 1 85%, equal 
to the calculated eclipse depth. 

Table 1. Relative photometry of WASP-4b. Extra digits are given 
to avoid round up errors. 



HJD 


Relative flux 


Uncertainty 


2455102.5150000 


0.99951 


0.00427 


2455102.5150069 


0.99732 


0.00427 


2455102.5150138 


1.00339 


0.00421 


2455102.5150277 


0.99452 


0.00424 


2455102.5150347 


1.00043 


0.00426 



Only a small portion of the data set is presented in this 
table, to exemplify its presentation format. The complete 
set can be found in the electronic version of the Journal. 

Table 2. Derived parameters of WASP-4b. 



Parameter 


Value 


68.3% Confidence Limits 


Unit 


Tc 


2455102.61162 


-0.00071, +0.00077 


HJD 


d 


0.185 


-0.014, +0.013 


% 


'^occ 


0.0881 


-0.0017, +0.0022 


d 


^coso; 


0.0027 


+ 0.0018 




Tb 


1995 


+ 40 


K 




Base-line parameters 




ao 


1.6933 


-0.0169, +0.0214 




ai 


-0.0364 


-0.0027, +0.0042 




a2 


-0.0112 


-0.0014, +0.0013 




as 


-0.00025 


-0.00074, +0.00062 




(34 


-0.00175 


-0.00060, +0.00078 





The center of the secondary eclipse occurs at = 0.49933 ± 
0.00059, assuming the ephemeris of Southworth et al. ( 20091) . 
For zero eccent ricity (Wilson et al. 2008; Gillon et al.l 120091: 
IWinn et al. 120091: ISouth worth et al 2009) . and a Hght travel time 
of 23.2 sec for this system (iLoebl 12005 ). the secondary eclipse 
is expected to occur at phase cpexp = 0.5002. The phase dif- 
ference S(p implies a non-zero eccentricity, that could be deter- 
mined knowing the prima ry transit length Ttr, from the equations 
(ICharbonneau et al.ll20Q5l) : 



ecosoj ^ n S(p , 



esmcj ^ 



Ttr + To, 



(3) 
(4) 
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The accuracy of our eclipse length measurement prevents us 
from using Eq. IH so we only use the former equation, to put a 
upper limit on the eccentricity e cos co - 0.0027 + 0.0018. Within 
the 3-(T level, this results argues in favor of a nearly circular 
orbit. 

Previous planet- to- star contrast measureme nts for WASP- 
4b have been obtained by iBeerer etal ] (120111) in the 3.6 and 
4.5 yum broad band filters with the Spitzer Space Telescope. 
As a zero-order approximation, we represent the planet with 
a black body w i th tem perature Tg, and the host star with a 
iHauschildt et aP (Il999l) stellar atmosphere model, with param- 
eters r,ff = 5500 /r, logg = 4.5, and [Fe/H] = 0.0 ( Gillon^etaLl 
l2009h . They are divided to obtain the planet- to- star flux con- 
trast as a function of wavelength. The expected depth was calcu- 
lated as weighted average of the contrast curve, weighted by the 
Ks filter transmission curve, the atmospheric transparency, and 
the detector efficiency curve. The best-fitting black body curve 
yielded a planetary brightness temperature of Tb - 1995 ± 40 K 
at ~ 2. 2 micron. 

The observational data of highly irradiated planets indicate 
very low Bond albedo value s (e.g. Charbonneau et al. 1999; 
iRowe et al.ll2008l: iRogers et aLll2009h . implying that the atmo- 
spheres of those planets are highly heated by the stellar ra- 
diation, and probably bloated. The temperature at which the 
planet re-emits the energy absorbed from the stellar flux, at a 
given distance a from the star, in units of the stellar radius, is 
given by Teq = T^^a^^^ (/(I - Ag))^/"^. Here, the / factor refers 
to the fraction of energy that is re-radiated from the planet to 
the observer. If the stellar radiation absorbed by the planet is 
re-radiated isotropically then / = 1/4, and we could deter- 
mine for a zero Bond albedo an equilibrium temperature for 
the planet of Teq = 1662 K. The difl'erence between the equi- 
librium and the brightness temperatures suggests a poor energy 
redistribution in the atmosphere of WASP-4b. In fact, to ob- 
tain the calculated brightness temperature the re-radiation fac- 
tor should be / = 0.52 , indicating a planet with re-radiation 
only from the day-side (iBurrows et al.l l2008l) . when assuming 
a zero Bond-albedo. T his result is in agreement with the results 
in lBeerer etaP (12011 ). who found that the Fortne v et all (l2008l) 
models with a small redistribution of heat, and an absence of 
TiO in the WASP-4b stratosphere, are best fits for their measure- 
ments. 

Planetary atmosphere models tuned for WASP-4b are shown 
in Fig. [21 The modeled atmospheres assume radiative and 
chemical equilibrium, and employ the chemical composi - 
tions and thermo-chemistry found in iBurrows &Sharpl (11999b . 
and Sharp & Burrows (2007), and the opacities described in 
[sharp & Burrows (2007), and references therein. These mod- 
els were calculated for diff'erent redistribution factors where 
Pn = implies no redistri bution, an d Pn = 0.5 means maximum 
redistribution (Burro ws et^ [20081 . This parameter represents 
the fraction of energy that is transferred from the day- side to the 
night- side of the planet. Moreover, some highly irradiated atmo- 
spheres require the inclusion of an extra absorber to fit calculated 
eclipse depths, this accounting for the presence of inversion lay- 
ers in the upper atmosphere ( Burro ws" et al.[[2007l [Knutson et al.[ 
[2QQ9[) . This scenario was tested by introducing an opacity pa- 
rameter Ke for this extra absorber as a second parameter to ex- 
plore. 

Our observations favor lower P^ parameter values, arguing 
in favor of an inefficient redistribution of the received radiation 
from the day- side to the night- side of the planet, as shown in 
Fig. [3 The models best suited with our measurement correspond 
to atmospheres with null or low amounts of a stratospheric ab- 



sorber, this leading to non inverted, or scarcely inv e rted a tmo- 
spheres. A similar result was found in "Beerer et alj ([201 1[) . for 
observations in the mid-IR. When considering the three mea- 
surements, i.e. the ISAAC one at 2.2 yum, and the Spitzer ones at 
3.6 and 4.5 yum, the best fitting models is the one with an ineffi- 
cient redistribution of heat, and an small quantity of absorber in 
the planet stratosphere. Thus, the lack of a strong thermal inver- 
sion in the atmosphere of WASP-4b could be inferred from this 
IR detection. 

4. Conclusions 

The ultrafast photometry technique developed in Paper I is able 
to detect very low amplitude secondary eclipses of extraso- 
lar planets. We have measured the secondary eclipse depth of 
WASP-4b which has an amplitude of = 0.185% at >10o- 
level. This secure detection strengthen our confidence in the new 
method, and indicates that it could be fruitful to attempt the ob- 
servation of secondary eclipses at shorter wavelengths, where 
the eclipses are shallower, but the NIR signal is still within our 
observational capabilities. 

For the stellar and planetary parameters, this yields a bright- 
ness temperature of T = 1995 K in the Ks, which agrees with 
the specific models for WASP-4b, and argues in favor of ineflft- 
cient heat redistribution from the day-side to the night-side of the 
planet. The absence of an strong thermal inversion in the strato- 
sphere of WASP-4b is inferred from our near-IR measurement, 
and mid-IR measurements in the literature. 
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Fig. 5. The planet- to- star flux contrast as a function of wavelength. The solid lines show three models with Kg = 0.00 - 0.03 cm^g, 
and Pn = OA - 0.3 for the planet WASP-4b. The squares are the estimated band- weighted measurements for these models, and the 
black point at 2.16yum represents our m easurement, with its error bars. Note that the red model point at 2.16 /urn is overlapped by 
our measurement. The measurements in iBeerer et al.lfcOl 1 ) are also shown at 3.6, and 4.5 jam. The transmission curves for the three 
filters are also shown for comparison as dotted lines. A color version of this figure is available in the electronic form. 
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